Violation of the Cauchy-Schwarz inequality with matter waves 
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The Cauchy-Schwarz (CS) inequality - one of the most widely used and important inequalities in mathematics 
- can be formulated as an upper bound to the strength of correlations between classically fluctuating quantities. 
Quantum-mechanical correlations can, however, exceed classical bounds. Here we realize four-wave mixing of 
atomic matter waves using colliding Bose-Einstein condensates, and demonstrate the violation of a multimode 
CS inequality for atom number correlations in opposite zones of the collision halo. The correlated atoms have 
large spatial separations and therefore open new opportunities for extending fundamental quantum-nonlocality 
tests to ensembles of massive particles. 
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The Cauchy-Schwarz (CS) inequality is ubiquitous in math- 
ematics and physics [1]. Its utility ranges from proofs of basic 
theorems in linear algebra to the derivation of the Heisenberg 
uncertainty principle. In its basic form, the CS inequality sim- 
ply states that the absolute value of the inner product of two 
vectors cannot be larger than the product of their lengths. In 
probability theory and classical physics the CS inequality can 
be applied to fluctuating quantities and states that the expecta- 
tion value of the cross-correlation (/1/2) between two quanti- 
ties Ii and I2 is bounded from above by the auto-correlations 
in each quantity: 
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This inequality is satisfied, for example, by two classical cur- 
rents emanating from a common source. 

In quantum mechanics, correlations can, however, be 
stronger than those allowed by the CS inequality [2-4]. Such 
correlations have been demonstrated in quantum optics using, 
for example, antibunched photons produced via spontaneous 
emission [5], or twin photon beams generated in a radiative 
cascade [6], parametric down conversion [7], and optical four- 
wave mixing [8]. Here flie discrete nature of the light and 
the strong correlation (or anticorrelation in antibunching) be- 
tween photons is responsible for the violation of the CS in- 
equality. The violation has even been demonstrated for two 
light beams detected as continuous variables [8]. 

In this work we demonstrate a violation of the CS inequal- 
ity in matter-wave optics using pair-correlated atoms formed 
in a collision of two Bose-Einstein condensates (BECs) of 
metastable helium [9-12] (see Fig. 1). The CS inequality 
which we study is a multimode inequality, involving integrated 
atomic densities, and therefore is different from the typical 
two-mode situation studied in quantum optics. Our results 
demonstrate the potential of atom optics experiments to ex- 
tend the fundamental tests of quantum mechanics to ensem- 
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Figure 1 . (Color online) Diagram of the collision geometry, (a) Two 
cigar-shaped condensates moving in opposite directions along the 
axial direction z shortly after their creation by a Bragg laser pulse 
(the anisotropy and spatial separation are not to scale), (b) Spheri- 
cal halo of scattered atoms produced by four-wave mixing after the 
cloud expands and the atoms fall to the detector 46 cm below. Dur- 
ing the flight to the detector, the unscattered condensates acquire a 
disk shape shown in white on the north and south poles of the halo. 
The (red) boxes 1 and 2 illustrate a pair of diametrically symmetric 
counting zones (integration volumes) for the average cross- and auto- 

f2) (2) 

correlation functions, Q12 and Qa (i — 1,2) (see text), used in the 
analysis of the Cauchy-Schwarz inequality. 



bles of massive particles. Indeed, violation of the CS inequal- 
ity implies the possibility of (but is not equivalent to) for- 
mation of quantum states that exhibit the Einstein-Podolsky- 
Rosen (EPR) correlations or violate a Bell's inequality [3]. 
The EPR and Bell-state correlations are of course of wider 
significance to foundational principles of quantum mechanics 
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than those that violate a CS inequality. Nevertheless, the im- 
portance of understanding the CS inequality in new physical 
regimes lies in the fact that: (;) they are the simplest possi- 
ble tests of stronger-than-classical correlations, and (//) they 
can be viewed as precursors, or necessary conditions, for the 
stricter tests of quantum mechanics. 

The atom-atom correlations resulting from the collision and 
violating the CS inequality are measured after long time- 
of-flight expansion using time- and position-resolved atom 
detection techniques unique to metastable atoms [13]. The 
307 ms long expansion time combined with a large collision 
and hence scattering velocity results in a ~6 cm spatial sep- 
aration between the scattered, correlated atoms. This sepa- 
ration is quite large compared to what has been achieved in 
recent related BEC experiments based on double-well or two- 
component systems [14—16], trap modulation techniques [17], 
or spin-changing interactions [18, 19]. This makes the BEC 
collisions ideally suited to quantum-nonlocality tests using ul- 
tracold atomic gases and the intrinsic interatomic interactions. 

In a simple two-mode quantum problem, described by bo- 
son creation and annihilation operators aj and a; (« = 1, 2), 
the Cauchy-Schwarz inequality of the form of Eq. (1) can be 
formulated in terms of the second-order correlation functions, 
G^^^ = (: fiihj :) = {a\a^-ajdi), and reads [2^] 

G^^ < [Gf'G^f (2) 

or simply G^^2 — ^ri symmetric case of G^^^ = ^^22- 

(2) (2) t 
Here, G\y = G^l^ hi = a\ hi is the particle number oper- 
ator, and the double colons indicate normal ordering of the 
creation and annihilation operators, which ensures the correct 
quantum-mechanical interpretation of the process of detec- 
tion of pairs of particles that contribute to the measurement 
of the second-order correlation function [2]. Stronger-than- 
classical correlation violating this inequality would require 
G^vl > [Gn G^2^]i/2, or g'^J > G^^^ in the symmetric case. 

The situation we analyze here is counterintuitive in that 
we observe a peak cross-correlation (for pairs of atoms scat- 
tered in opposite directions) that is smaller than the peak auto- 
correlation (for pairs of atoms propagating in the same direc- 
tion). In a simple two-mode model such a ratio of the cross- 
correlation and auto-correlation satisfies the classical CS in- 
equality. However, in order to adequately treat the atom-atom 
correlations in the BEC collision problem, one must general- 
ize the CS inequality to a multimode situation, which takes 
into account the fact that the cross- and auto-correlations in 
matter- wave optics are usually functions (in our case of mo- 
mentum). The various correlation functions can have differ- 
ent widths and peak heights, and one must define an appro- 
priate integration domain over multiple momentum modes to 
recover an inequality that plays the same role as that in the 
two-mode case and is actually violated, as we show below. 

The experimental setup was described in Refs. [11, 12]. 
Briefly, a cigar-shaped BEC of metastable helium, contain- 
ing approximately ^ 10^ atoms, trapped initially in a har- 
monic trapping potential with frequencies (a;^., w^, aj2)/27r — 



(1500, 1500, 7.5) Hz, was split by Bragg diffraction into two 
parts along the axial (z-) direction [see Fig. 1(a)], with ve- 
locities differing by twice the single photon recoil velocity 
Vice — 9.2 cm/s. Atoms interact via binary, momentum con- 
serving s-wave collisions and scatter onto a nearly spherical 
halo [see Fig. 1(b)] whose radius in velocity space is about the 
recoil velocity [11, 20]. The scattered atoms fall onto a detec- 
tor that records the arrival times and positions of individual 
atoms [13] with a quantum efficiency of ~ 10%. The halo di- 
ameter in position space at the detector is ^ 6 cm. We use the 
arrival times and positions to reconstruct 3D velocity vectors 
V for each atom. The unscattered BECs locally saturate the 
detector To quantify the strength of correlations correspond- 
ing only to spontaneously scattered atoms, we exclude from 
the analysis the data points containing the BECs and their im- 
mediate vicinity (|u^| < 0.5 Wroc) and further restrict ourselves 
to a spherical shell of radial thickness 0.9 < Vr/vi-cc < 11 
(where the signal to noise is large enough), defining the total 
volume of the analyzed region as Vdata- 

Using the atom arrival and position data, we can mea- 
sure the second-order correlation functions between the atom 
number densities f;,(k) at two points in momentum space, 
(7(^)(k, k') = (: n(k)7i(k') :) (see Supplementary Material 
[21]), with k denoting the wave-vector k ~ mv/fi, and tik the 
momentum. The correlation measurements are averaged over 
a certain counting zone (integration volume V) on the scatter- 
ing sphere in order to get statistically significant results. By 
choosing k' to be nearly opposite or nearly collinear to k, we 
can define the averaged back-to-back (BB) or collinear (CL) 
correlation functions, 

g^iiAk)^ [ d'kg^'\k,~k + Ak), (3) 

4'i^(Ak)= / d'kg^^Hk,k + Ak), (4) 

JV 

which play a role analogous to the cross- and auto-correlation 

(2) (2) 

functions, G\2 ™d G^ , in the simple two-mode problem 
discussed above. The BB and CL correlations are defined as 
functions of the relative displacement Ak, while the depen- 
dence on k is lost due to the averaging. 

The normalized BB and CL correlations functions, 
.gg^g(Ak) and g'^^{Ak), averaged over the unexcised part of 
the scattering sphere Vdata are shown in Fig. 2. The BB cor- 
relation peak results from binary, elastic collisions between 
atoms, whereas the CL correlation peak is a variant of the 
Hanbury Brown and Twiss effect [22, 23] — a two-particle 
interference involving members of two different atom pairs 
[9, 10, 24, 25]. Both correlation functions are anisotropic be- 
cause of the anisotropy of the initial colliding condensates. 

An important difference with the experiment of Ref. [9] is 
that the geometry in the present experiment (with vertically 
elongated condensates) is such that the observed widths of the 
correlation functions are not limited by the detectorresolution. 
Here we now observe that the BB and CL correlations have 
very different widths, with the BB width being significantly 
larger than the CL width. This broadening is largely due to the 
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Figure 2. (Color online) Normalized back-to-back (a) and collinear 
(b) correlation functions, (?g^g(Ak) and gQl{A\i), in momentum 
space integrated over Vdata coiTesponding to \kz\ < 0.5 fcroc and 
0.9 < kr/kicc < 1.1, where k rcc — TTlVrcc /h is the recoil momen- 
tum. The data is averaged over 3600 experimental runs. Because of 
the cylindrical symmetry of the initial condensate and of the overall 
geometry of the collision, the dependence on the k^c and ky compo- 
nents should physically be identical and therefore can be combined 
(averaged); the correlation functions can then be presented as 2D sur- 
face plots on the (kz, k^y) plane. The 2D plots were smoothed with 
a nearest neighbor running average. The data points along the kz and 
kxy projections (coiTesponding to thin slices centered at k^y ~ 
and kz = 0, respectively) are not smoothed. The solid lines show 
the Gaussian fits to these projections. The peak height of the back- 
to-back correlation function is ~ 1.2 while that of the collinear corre- 
lation function is ~ 1.4, apparently confirming the Cauchy-Schwarz 
inequality. The widths of the two distributions are, however, very 
different (aBB,x — o"bb,j; — 0.21fcrcc, ctbb.z — 0.019/crcc, whereas 
acL,x — o'CL.y — 0.036fcrec, crcL,2 — 0.002fcrcc) and a multimode 
formulation of the Cauchy-Schwarz inequality, which relates the rel- 
ative volumes of the correlation functions, is violated. 

size of the condensate in the vertical direction (^ 1 mm). The 
elongated nature of the cloud and the estimated temperature 
of ^ 200 nK also means that the condensates correspond in 
fact to quasicondensates [26] whose phase coherence length 
is smaller than the size of the atomic cloud. The broadening 
of the BB correlation due to the presence of quasicondensates 
will be discussed in another paper [27], but we emphasize that 
the CS inequality analyzed here is insensitive to the detailed 
broadening mechanism as it relies on integrals over correla- 
tion functions. This is one of the key points in considering the 



multimode CS inequality. 

Since the peak of the CL correlation function corresponds 
to a situation in which the two atoms follow the same path, 
we can associate it with the auto-correlation of the momen- 
tum of the particles on the collision sphere. Similarly, the 
peak of the BB correlation function corresponds to two atoms 
following two distinct paths and therefore can be associated 
with the cross-correlation function between the respective mo- 
menta. Hence we realize a situation in which one is tempted 
to apply the CS inequality to the peak values of these corre- 
lation functions. As we see from Fig. 2, if one naively uses 
only the peak heights, the CS inequality is not violated since 
A^) < g^cliO) and hence ^(f^(0) < gg^(O) due to the 
nearly identical normalization factors [21]. 

We can, however, construct a CS inequality that is violated 

(2) 

if we use integrated correlation functions, Q^j , that corre- 
spond to atom numbers Ni = (PIl a'' (k) a(k) {i ~ 1, 2) 
in two distinct zones on the collision halo [21]: 

= (: N,N, :)= [ d'k [ d'k' G^^HK^')- (5) 
Jvi Jv, 

The choice of the two integration (zone) volumes Vi and Vj 

_(2) 

determines whether the Q^j -function corresponds to the BB 
{i 7^ j) or CL (z = i) correlation functions, Eqs. (3) and (4). 
The CS inequality that we can now analyze for violation 

reads g^^^ < {QllQ-^^f^'^ ■ To quantify the de gree of viola- 
tion, we introduce a correlation coefficient 

-(2)-(2) 1/2 
C = ^22 J , (6) 

which is smaller than unity classically, but can be larger than 
unity for states with stronger-than-classical correlations. 

In Fig. 3 we plot the correlation coefficient C determined 
from the data for different integration zones Vi and V2, but al- 
ways keeping the two volumes equal. When Vi and V2 corre- 
spond to diametrically opposed, correlated pairs of zones (red 
circles), C is greater than unity, violating the CS inequality, 
while for neighboring, uncorrected pairs (blue squares) the 
CS inequality is not violated. The figure also shows the results 
of a quantum-mechanical calculation of C using a stochastic 
Bogoliubov approach (green solid curve) [20, 21, 28]. The 
calculation is for the initial total number of atoms iV = 85 000 
and is in good agreement with the observations. The choice of 
large integration volumes (small number of zones M) results 
in only weak violations, while using smaller volumes (large 
M) increases the violation. This behavior is to be expected 
[21] because large integration zones include many, uncoiTe- 
lated events which dilute the computed correlation. The sat- 
uration of C, in the current arrangement of integration zones 
- with a fixed number of polar cuts and hence a fixed zone 
size along z which always remains larger than the longitu- 
dinal correlation width - occurs when the tangential size of 
the zone begins to approach the transverse width of the CL 
correlation function. If the zone sizes were made smaller in 
all directions, we would recover the situation applicable to the 
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Figure 3. (Color online) Correlation coefficient C as a function of 
the number of zones M = Vdata/Vi into which we cut the scatter- 
ing sphere. C > 1 corresponds to violation of the Cauchy-Schwarz 
inequality. The scattering sphere was cut into 8 polar and from 2 to 
80 azimuthal zones; the resulting arrangement of zones for M — 16 
and 32 is illustrated in the upper panel. The observed values of C 
for pairs of correlated diametrically opposite zones (shown in red in 
the upper panel as an example) were averaged to get one data point 
for a given M; the data points for such zones are shown as red cir- 
cles, for uncoiTelated (neighbouring) zones — as blue squares. The 
error bars show the standard deviation of the mean over the number 
of zone pairs. The (green) solid curve is the theoretical prediction 
[21] calculated using the experimental parameters and a stochastic 
Bogoliubov approach [20, 28]. 

peak values of the correlation functions (and hence no CS vio- 
lation) as soon as the sizes become smaller than the respective 
correlations widths (see Eq. (Sll) in [21]). 

We have shown the violation of the CS inequality using the 
experimental data of Ref . [ 1 1 ] in which a sub-Poissonian vari- 
ance in the atom number difference between opposite zones 
was observed. Although the two effects are linked mathe- 
matically in simple cases, they are not equivalent in general 
[8, 21]. Because of the multimode nature of the four-wave 
mixing process, we observe stronger (weaker) suppression of 
the variance below the shot-noise level for the larger (smaller) 
zones (see Fig. 3 of [11]), whereas the degree of violation of 
the CS inequality follows the opposite trend. This difference 
can be of importance for other experimental tests of stronger- 
than-classical correlations in inherently multimode situations 
in matter-wave optics. 

The nonclassical character of the observed correlations im- 
plies that the scattered atoms cannot be described by classi- 
cal stochastic random variables [29]. Our experiment is an 
important step towards the demonstrations of increasingly re- 
strictive types of nonlocal quantum correlations with matter 
waves, which we hope will one day culminate in the viola- 
tion of a Bell inequality as well. In this case, the nonclassi- 
cal character of correlations will also defy a description via 
a local hidden variable theory [4, 29]. Non-optical violations 
of Bell's inequalities have so far only been demonstrated for 
pairs of massive particles (such as two trapped ions [30] or 
proton-proton pairs in the decay of ^He [31]), but never in the 



multi-particle regime. The BEC collision scheme used here is 
particularly well-suited for demonstrating a Bell inequality vi- 
olation [32] using an atom optics analog of the Rarity-Tapster 
setup [33]. 
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